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Introduction 

NmerouJ landfills and Industrial activities In North America are 
situated such that chemicals emanating from these facilities are 
contaminating groundwater in underlying bedrock. In particular, this 
Is the ease for sites in southern Ontario and the Niagara Peninsula 
where the underlying materials consist of a thin veneer of overburden 
overtop moderately to highly fractured shale of the Ordovlclan 
Queenston Formation or the Cambrian Georgian Bay Formation. Because 
contaminant migration along fracture planei within the shale can be 
very rapid relative to unconsolidated media, accurate prediction of 
the direction, velocity and dilution of the contamination Is 
especially Important with regard to public health and safety. 

Groundwater velocity predictions are usually based on direct 
hydraulic measurements of the permeability of fractures which is 
expressed as an equivalent single fracture aperture, 2b. Average 
groundwater velocity, v. Is then estimated from the expression: 



12)1 di 



m 



where p. g and p are the density, gravitational acceleration and 
dynamic viscosity of the water, respectively. The hydraulic gradient, 
dh/dl. Is usually deterained frcoi water level measurements tn plezo- 
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»eurs or more correctly from w)tf-i,vel c.sing string, which IsoUt. 
individual fractures and fracture zones. 

Equation (I) is derived from a conceptual model known as the 
cubic lav 1n which water flow ,s between two sn»oth parallel plate, 
separated by the width. 2h. Inspection of natural fractures In 
outcrop and core sample show that fracture surfaces are distinctly 
non-uniform and are characterized hy a variety of aperture widths 
regions where asperities contact and channels of reduced or enhanced 
opening. Consequently, the relationship between the 2b obtained from 
the results of hydraulic tests Interpreted using the cubic law and the 
-erage mechanical or true aperature Is unclear. This Is evidenced by 
ccmpanng the groundwater velocity predicted fr«. apertures obtained 

ZTc, T"'!" ™""^'""' '" ""' »"»'-'' f™ tracer experi™„ts 
wMch Show that hydrauHcally measured apertures (hereafter called 

,ZT\r"""' "'^"*'-'« groundwater velocity In .ost cases 
(Novakowski et al.. . isaS; Raven et al-. 1988). 

The Objective of this study ls to investigate the relation 
between fracture apertures determined from hydraulic methods versus 
those determined using tracer experiments with the Intent to resolve 

If". ' """"'' "'" " '''" ^'""^ '"'""ing on a 

ins e-flat lying fracture Intersected by four boreholes of a seven 

steady-state single-well and multiple-well hydraulic tests were 
employed to determine the hydraulic aperture. The tracer apertures 
-ere determined from seven tracer experiments conducted using both 
mjectlon-wuhdrawal and radl.l-convergent flow fields configurations. 

Field Method 

whlcht'^! T ''''"' °' '''' ""^ ''''■ ' ''*"" =''^ -" =^»"n at 
Which preliminary hydraulic tests revealed the presence of two largely 
unconnected fl,t-.ylng fractures at depths of 9.7 and 10.4 m 
groun surface. The study site Is located ,n Clarkson, Ontario 
unus. petro-canada property about l.S km north of the Lake Ontar 
.horellne. The fractures are In the Meaford-Oundas shale whic 



274 



sub-crops In the Oakvllls-Hlsslssauga area of southern Ontario and 
strikes North-Northwest. 

The Initial location of the fracture set and the subsequent 
characterization of the apertures and boundaries were carried out over 
198G and 1967. During 19S6. a preliminary drilling program was 
conducted In which five 76 mm boreholes were drilled IS m apart 1n two 
perpendicular rows of three (Figure 1). Examination of the core and 
preliminary constant-head test results showed four main fracture zones 
with the shallowest set at 9.5 m to 10.5 m exhibiting the highest 
permeability. Subsequent field work during 1987 focused on this set 
and two new boreholes, UW6 and UH7 were drilled during mid-sunner 
solely to Investigate this zone. After completion of the drilling. 
the precise location of each of the two fractures in all seven bore- 
holes was deteniilned using constant-head tests with a O.I m Interval 
spacing. 

On the basis of the constant-head test results, the upper frac- 
ture has a mean aperture width of about 120 pm (this is equal to an 
hydraulic conductivity of 1 x lO'^ m/s) and the lower fracture has a 
mean aperture of about 200 \m (hydraulic conductivity of 
5 X Ifl-S m/s). Figure 1 shows the approximate extent of each fracture 
wUhln the study site. The boundaries Indicated on Figure 1 are arbi- 
trary and are drawn only to show Interconnection between boreholes. 
For example, fracture no. 1. the upper fracture, hydraullcally 
connects UW3, UH6 and UW7 but Is not well connected to UWJ. UW2, UW4 
and UW5. On the basis of the larger aperture Identified for the lower 
fracture, it was selected for the purpose of this study and all 
further hydraulic and tracer tests described 1n the following text 
were conducted In this fracture. Preliminary estimates of hydraulic 
head In the lower fracture indicate groundwater flow oriented towards 
the northeast with a gradient of 0.01. 

In addition to constant-head tests, slug tests and steady-state 
and transient pumping tests were also conducted. The slug tests were 
carried out using a straddle-packer arrangement attached to a 1.5" 
I.D. standpipe in which the Instantaneously Induced rise or fall of 
water level was »on1tored using a pressure transducer. Slug tests 
were conducted in all boreholes except UW6. The pumping tests were 
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FIELD STU3Y SITE 
Ordovician Shale 




FractuTB De oths: 
Na 1 - 9.7 mbfls 
Na 2 - 10.4 mbgs 

Na 1 > 75 iim 

No. 2 > 175 nm 



Figure 1. Approximate boundaries of the two fractures Intersected by 
boreholes at the study site. Fracture no. 1 Interconnects 
boreholes UM3, UW6 and UW7. Fracture no. Z Interconnects 
UWl, UU2, UW4, and UW5. 
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conducted also using a straddle packer arrangenwnt except in this esse 
the test fracture was completely Isolated. Abstraction of water was 
carried out by pumping from surface through f nylon tubing. Drawdown 
In the Observation Intervals also Isolated by packers was monitored 
using pressure transducers. Pumping tests were conducted In UW5 and 
UW2 with observation In UWl, UW2 and iUV4. 

Two methods for Induced-gradient tracer experiments were en^loyed 
to compare to the resolts of the hydraulic tests; the Injection- 
withdrawal loethod and the radial-convergent ™thod. The injection- 
withdrawal tracer experiments were conducted by Injecting water Into 
one borehole and withdrawing water from another at the sa»e flowrate. 
Sodium bromide tracer was Introduced as a slug Into the Injection 
borehole once the flow field had achieved steady-state conditions. 
Tracer arrival at the withdrawal borehole was continuously monitored 
using a flow-cell arrangeaent and tracer concentration was measured In 
discrete samples of the withdrawn water. The fracture was Isolated In 
both the Injection and withdrawal boreholes using straddle packer 
arrangements. Injection and withdrawal flowrates and Interwell 
hydraulic gradients were monitored on a continuing basis through the 
duration of each experiment. The radial-convergent tracer experiments 
were conducted In a similar manner except that continuous Injection of 
water was not .employed. A steady radial flow field was established 
around a borehole by pusupmg and tracer was Introduced Into the flow 
field from an observation borehole by pulse Injection. Both tracer 
methods were employed for tests conducted between UWl and UW2 and UW4 
and UMZ. with UW2 as the withdrawal borehole in all cases. The 
mjectlon-wlthdrawal method only was employed for experiments 
conducted between UWZ and UW5, UM5 and UW4 and UW7 and UW3 
NovakowskI (1988) provides a more detailed description of the field 
methods employed for this study. 

Results and Discussion 

The constant-head hydraulic tests were analysed using the Thiem 
equation where the steady flow rate. 0, divided by the rise In 
hydraulic head, iH, Is proportional to transmtsslvlty. T: 
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(2) 



where r^ Is the radtus of tnfluenM (estimated at approximately 10 m 
for these tests) and r« Is the radius of the well. Fracture 
aperture Is then calculated from the trarsralsslvity using: 



2b 



fiZU ,,1/3 



(3) 



Both the slug tests and transient pumping tests were analysed 
using type curve methods to obtain transralsslvlty. The steady-state 
pun.plng tests were analysed using equation (2) where the radial 
disunce between the pumping and observation boreholes Is substituted 
for re and 4H Is the difference In hydraulic head between the two 
wells. 

The tracers experiments were analysed for fracture aperture by 
using analytical approximations to describe the streM.1 1ne(s) between 
the injection and pumping borehole and solving the advectfon-dlsper- 
slon equation along each streamline. The analytical approximations 
for the streamlines are based on the average travel time for the 
tracer and the interwell difference In hydraulic head. The aperture 
is determined by a trial and error procedure, varylns the Zb to 
determine the location of the model peak. 

Table I shows the results of the hydraulic tests and tracer 
experiments expressed as fracture aperture In micrometers along with 
the borehole Interval Isolated for each test. In general, distinction 
between the results of transient and steady-state hydraulic tests is 
immediately clear. The results of the transient hydraulic tests show 
consistently larger fracture apertures than those of the steady-state 
tests. This is a function of the radius of test Influence and knowing 
that steady-state test results more likely reflect the properties of 
the fracture nearest the source borehole, the apertures deteralned 
from the transient tests probably reflect regions of the study 
fracture that have a slightly larger aperture and that are beyond the 
lateral extent of the study site. 



278 



?_iriL- 




On the basis of the relative unlfonnUy of the apertures deter- 
mined from the transient hydraulic tests conducted In borehole WC 
{the central borehole of the study site) there Is no evidence to 
suggest non-radial flow effects which might Indicate the presence of 
channeled flow. This suggests that the points of closure currently 
Identified by boreholes UW3. UH6 and U¥7 are relatively small In area 
and the boundaries to the fracture plane as drawn In Figure 1 are 
probably inappropriate. 

Comparing the fracture apertures determined from the hydraulic 
tests to those from the tracer experiments. It Is evident that the 
tracer apertures are consistently smaller by up to SOS. This observa- 
tion Is In agreement with the results of other field experiments 
reviewed In Raven et al. (1988). These results are significant with 
respect to predicting groundwater velocity In fractures because 
velocity Is proportional to the square of the fracture aperture. 
Therefore, an aperture difference of a factor of 2, for example, 
translates Into a four-fold difference In velocity. 

An explanation as to why tracer apertures are always snaller than 
hydraulic apertures can be found based on the solute transit path 
length. Because we have established that natural fractures are rough 
and undulatlngi, and we know that the fracture walls are in some degree 
of contact as evidenced by closure In UW3, l(H6 and UW7, then the path 
travelled by the tracer In transit from the Injection borehole to the 
pumping borehole may be very tortuous as the flow negotiates around 
the points of closure and areas of reduced opening. 

To investigate this explanation, the tracer experiments between 
UWl and UH2 and UW4 and LIW2 were simulated by fixing the 2b to the 
size of the steady-state hydraulic aperture and modeling for the path 
length. Table 2 shows the modeled path lengths and values of 
tortuosity which are calculated by dividing the modeled path length 
Into the actual Interwell distance. 
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Table I. Aperture widths (In pm) as determined from slufl test, 
constant head test, tracer test and pumping test results. 



Borehole 



51 ua Tests 

UWl 
UKZ 

im 

UMI: 

MS 
M? 

Constant Head Tests 

UWl 
UM2 
IM3 

UW4 
UW5 
UW6 
UW7 

Tracer Test? 



UM1-UW2(I 

UW1-UW2(P 

UH2-UW5(I 

UW4-UW2(Ij 

UW4-UW2(I>) 

UV(5-*UW4(l5 

uw;-uw3(i) 



Pumping Tests 

UW5-UW2(T) 
UW5-UWl(Ti 
UWZ-UW4(S) 
UW2-UW1(S) 
UW2-rtJWl(T) 
UM2-UV(4(T) 



/^trture 



210 

267 
147 
226 

180 

140 



229 
236 
13S 

1S6 
164 
80 
120 



160 
132 
192 
1S5 
ISI 
164 
151 



184 
179 
211 
215 
254 
228 



T - transient 

S - Steady state 

I - tnjectlon-wlthdramal flow field 

P - radial-convergent flow field 



Interval 
(mbgs) 



10.10 
10.16 
9.40 
10.21 
10.20 



10.92 
10. 9B 
10.22 
11.03 
11.02 



9.90 - 10.52 



9.90 - 


10.72 






9.96 - 


10.76 






9.10 - 


9.92 






10.21 - 


11.03 






9.80 - 


10.62 






9.60 - 


10.42 






9.30 - 


10.12 






10.10 - 


10.92 


10.00 - 


10.49 


10.10 - 


10.92 


10.00 - 


10.49 


10.00 - 


10.49 


10.00 - 


10.90 


10.20 - 


11.02 


10.00 - 


11.49 


10.20 - 


11.02 


10.00 - 


10.49 


10.40 - 


10.89 


10. ZO - 


11.02 


9.50 - 


9.99 


9.40 - 


10.22 


9.80 - 


11.02 


10.10 - 


10.95 


9.80 - 


11.02 


10.10 - 


10.95 


10.00 - 


10.49 


10.10 - 


10.92 


10.00 - 


10.49 


10.20 - 


11.02 


10.00 . 


10.49 


10.10 - 


10.92 


10.00 - 


10.49 


10.10 - 


10.92 



r 



Table 2. Modeled 

tracer 


path lengths and tortuosity 
experiments. 


determined from the 






Path Length (m) 


Tortuoustty 


UHl to UW2(I) 
UWl to UW2(P) 
UW4 to UWZ(I) 
UW4 to UW2(P) 




19.0 
23.0 
20.6 
20.6 


0.79 
0.65 
0,73 
0.73 



The range of tortuosity Is from 0.65 to 0.79 which is quite wide 
although not unrealistic In terms of the average angularity of flow 
(the vectoraV direction of flow not directly towards the withdrawal 
well). Tortuosity of 0.65 gives an angularity of flow of 50* while 
D.79 gives about 38'. The angle of flow In porous media Is about 
45". Considering the complex possibility of pathways dependent on the 
aperture distribution, these vtlutt of tortuosity are probably 
reasonable. 

On the basis of the results presented here and on other recent 
field evidence (Raven et al. . 1988), It Is apparent that using the 
cubic law to interpret the results of hydraulic tests will lead to 
erroneous predictions of groundwater velocity. Therefore a new 
conceptual model Is urgently required with which the results of 
hydraulic tests can be utilized to predict velocity without relying on 
coiBparlson to more costly and time consuming tracer experiments. On 
the basis of the preliminary modeling conducted for this study. It Is 
tikely that such a conceptual model will depend on the role of 
tortuousity during solute transport In natural fractures. 

Sumiary 

Comparisons of fracture aperture widths determined from the 
results of hydraulic tests and tracer experiments were made. 
Constant-head, slug, steady-state and transient pwnping tests were 
employed to determine the fracture aperture, hydraullcally. The 
tracer experiments were carried out In radial-convergent and 
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Injection-Withdrawal flow field formats. The tests were conducted In 
a discrete flat-lying fracture In Ordovlclan aged shale at « depth of 
approximately 10 m. 

The results show that apertures determined from tracer 
experiments are consistently smaller than those obtained from 
hydraulic tests. The reason for this discrepancy Is probably due to 
the Influence of tortuous flow around reduced apertures and points of 
closure In the fracture plane. The development of a new conceptual 
model which can more accurately relate the results of hydraulic tests 
to a prediction of groundwater velocity, ts required to reconcile 
these differences. 
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